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Abstract 

We present results for the one loop anomalous dimension matrix of flavor non-singlet twist-4 
operators of lowest spin that contribute to the leading moment of the F2 structure function in deep 
q_j inelastic electron-nucleon scattering. We analyze the flavor structure of the anomalous dimension 

matrix and decompose the leading moment of F2 into separate flavor channels. In addition to 
building on previous work with higher-twist operators, these results can provide a benchmark for 
future work that generalizes to include the higher moments as well. We include non-perturbative 
input from the lattice and phenomenological estimates of the twist-4 matrix elements and estimate 
the twist-4 contributions to the leading moment of Fi- The results suggest that the overall twist-4 
contribution may be suppressed due to either cancellations among the twist-4 terms or inherently 
small twist-4 matrix elements. 
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I. INTRODUCTION 



In this paper we report on a computation of the leading log Q 2 evolution of the twist- 
4 contribution to the flavor non-singlet, leading moment of the F 2 structure function and 
provide phenomenological estimates of the resulting effects. Our study is motivated broadly 
by one of the major challenges for nuclear physics: understanding the dynamics of quarks 
and gluons, as determined by the QCD Lagrangian, and explaining their connection to 
the hadronic degrees of freedom. The twist expansion in QCD is a useful tool that is 
well-suited for this challenge. The property of asymptotic freedom of QCD allows one to 
calculate sufficiently inclusive hadronic observables, at asymptotically high energies, in terms 
of the perturbative quark and gluon degrees of freedom. This phenomenon is often referred 
to as quark-hadron duality. One of the simplest examples of this duality is the process 
e + e~ — > hadrons which is described well by the quark-level process e + e~ —> qq away from 
thresholds. Similarly, deep inelastic electron-nucleon scattering is described by electrons 
scattering off free quarks in the asymptotic region. 

However at low momentum scales, where the strong coupling a s is large, multi-parton 
correlations become important and lead to violations of quark-hadron duality. These correla- 
tions are embodied in higher twist effects. At sufficiently low scales, QCD is non-perturbative 
and the hadronic bound states of quarks and gluons become the relevant degrees of freedom. 
Despite this, there are several examples where low energy hadronic observables, averaged 
over appropriate intervals, exhibit behavior that reveal the underlying connection to the 
quark and gluon degrees of freedom. Bloom and Gilman [1-3] first observed that electron- 
nucleon scattering in the resonance region is related to the deep inelastic scaling regime. 
In particular, they observed that even in the region of low momentum transfer (Q 2 ), domi- 
nated by the highly non-perturbative dynamics of nucleon resonances, the nucleon structure 
function F 2 (x,Q 2 ) exhibits logarithmic scaling in Q 2 when averaged over appropriate inter- 
vals in Bjorken-x. Furthermore, the resonance structures seen in F 2 (x,Q 2 ) as a function 
of the Bjorken-x, slide along the deep inelastic scaling curve for increasing Q 2 . This loga- 
rithmic scaling in Q 2 of the structure function F 2 is described by the DGLAP evolution of 
the leading twist parton distribution functions (PDFs). This manifestation of quark-hadron 
duality which relates the resonance region to the deep inelastic scaling region is known as 
the Bloom-Gilman (BG) duality. At low values of Q 2 , one expects the onset of power law 
behavior corresponding to contributions from higher twist terms in the operator product 
expansion (OPE). Such behavior would signal a clear violation of BG duality and give a 
direct probe of multi-parton correlations in the nucleon. 

Detailed studies of the BG duality and its violation can provide insight into the dynamics 
of the quark-hadron transition and have lead to a large experimental effort. Since the early 
days of the SLAC-MIT [3] experiment, a wealth of data (for a comprehensive review see 
Ref.[4]) on structure functions has been accumulated over a wide range in x and Q 2 . A 
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large fraction of this data [5-11] is on the proton structure function F 2 (x,Q 2 ) which is 
now the best measured quantity in deep inelastic electron scattering. Data is also available 
on deuterium [12] and heavy nuclear [13] targets in the high-x and low Q 2 region. The 
data on the structure functions in the resonance region has been compared to the scaling 
curves obtained from global fits [14, 15] of the PDFs and DGLAP evolution. The resonance 
peak structures seen in the F 2 (x, Q 2 ) structure function are observed on average to oscillate 
around the scaling curve. In particular, the average of the structure function over all values 
of x, including over all resonance peaks, exhibits scaling behavior. Furthermore, the average 
of F$(x, Q 2 ) over individual resonance peaks is also observed to follow the scaling curve and 
is known as local BG duality. Similar scaling behavior is also observed for deuterium and 
heavy nuclei targets. 

In modern field-theoretic language, quark-hadron duality can be quantitatively formu- 
lated [16, 17] in terms of the operator product expansion (OPE). The Cornwall- Norton 
moments of the F 2 nucleon structure function 

Mt\Q 2 )= f dxx n - 2 F 2 (x,Q 2 ), (1) 
J o 

can be expressed schematically in terms of the OPE as 

Mt\Q 2 ) = E E (^V C nrM 2 ) ( 2 ) 
i r=2k W / 

Here k runs over all positive integers and the indices n, r, i denote the spin, twist, and type 
of operators respectively. The twist is defined r = d — s, where d, s denote the dimension and 
spin of the operator O l nT . The Wilson coefficients C l nT are perturbatively calculable as an 
expansion in a s (Q 2 ) and exhibit logarithmic scaling in Q 2 /fJ, 2 with /j, being an appropriately 
chosen input scale. The non-perturbative nucleon matrix element of the operator O l nT is 
denoted by (O l nr ) and has been scaled to an appropriate power of a typical hadronic scale 
A ~ 1 GeV. The power law behavior in Q 2 of the various terms in the OPE is determined by 
the twist r. The leading twist (r = 2) nucleon matrix elements are given by the moments of 
the standard PDFs. Quark-hadron duality corresponds to the dominance of the leading twist 
terms which are determined by the scattering of electrons from almost free quarks weighted 
by the PDFs. Logarithmic corrections to Bj or ken-scaling are determined by the standard 
DGLAP evolution of the PDFs. Violations of quark-hadron duality arise from the higher 
twist terms in the OPE as power corrections in 1/Q 2 . These higher twist terms encode long 
range multi-parton correlations in the nucleon and are expected to become important at low 

Q 2 . 

In the language of the OPE, the observed BG duality corresponds to unexpectedly small 
contributions from the higher twist terms to the lowest (n = 2) moment of the F 2 struc- 
ture function at low Q 2 . The higher moments of the structure function, weighted more by 
the large x resonance region, are expected to be more sensitive to higher twist effects. A 
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recent analysis by the CLAS collaboration [11] found that the moments of the F%(x, Q 2 ) 
structure function were dominated by the leading twist terms down to Q 2 ~ 1 GeV 2 , im- 
plying correspondingly small higher twist effects. For the lowest moment, after accounting 
for kinematic power corrections, the higher twist contributions were less than about 5% of 
the leading twist moment for Q 2 > 1 GeV 2 . These results were obtained through a detailed 
study of the Q 2 behavior of the collected data. 

The high quality of available data allows for a systematic study of higher twist correla- 
tions, providing a window into quark-hadron duality violations and nucleon structure. One 
limitation for such a program is the lack of precise theoretical knowledge of the renormaliza- 
tion group (RG) evolution of the higher twist operators. Given the absence of this theoretical 
input, the CLAS collaboration considered the effects of twist-4 and twist-6 contributions, in 
addition to the leading twist effects, using a simple ansatz [18, 19] that parameterizes these 
contributions to the moments of the structure function as 



Mt\Q 2 )=Vn(Q 2 ) + a^ 



a s {Q'' 



,2 



^- + a (6) 



a s {Q 2 



(6) 

2\"|7n ,,4 



where r] n (Q 2 ) is the leading twist contribution, and parameterize the twist-4 and 
twist-6 nucleon matrix elements respectively, and "fn^ and are effective anomalous dimen- 
sions parametrizing the RG evolution of twist-4 and twist-6 operators respectively. With the 
parameterization written in Eq.(3), the CLAS collaboration interpreted the unexpectedly 
tiny higher twist contribution to the moment as being due to a conspiracy of cancellation 
between twist-4 and an oppositely signed twist-6 contribution. From a rigorous theoretical 
perspective, however, the situation is considerably more complex, as higher twist contribu- 
tions are determined by a large number of operators that mix under RG evolution. Both 
the contributions from these matrix elements and the details of their mixing are ignored in 
Eq.(3). A basis of operators at twist-4 along with their tree level Wilson coefficients was 
first given in Refs. [20-22], in the transverse basis in Ref. [23], and more recently using the 
soft-collinear effective theory (SCET) in Ref. [24]. A conformal basis of higher twist opera- 
tors was constructed in Ref. [25] and a one-loop analysis of conformal higher twist operators 
is presented in Refs. [25, 26]. 

Higher twist operators in QCD are also of interest for parity violating deep inelastic 
scattering (PVDIS). As part of the 12 GeV upgrade at JLAB, new experiments [27, 28] will 
measure the electron polarization asymmetry, 

+ Oh 

in parity violating deep inelastic scattering off a deuteron target over a wide range of Q 2 
and x to sub-precent level precision. The impact of the logarithmic running of higher twist 
operators is one effect one must account for when interpreting the asymmetry. Due to the 
high precision of the measurements, hadronic uncertainties including higher-twist effects 
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must be investigated carefully as they can potentially cloud theoretical interpretations of 
deviations from Standard Model (SM) predictions. The effects of higher twist contributions 
to this parity- violating asymmetry were recently studied in Refs. [29, 30]. Based on the 
argument by Bjorken [31] and Wolfenstein [32] it was shown [30] that this asymmetry can 
be a powerful probe of quark-quark correlations in the nucleon. For a deuterium target, Arl 
is sensitive to a single four-quark operator involving up and down-quark fields 



which is a twist-4 operator. As pointed out in [30], combining high precision data taken 
over a wide range of x and Q 2 from future PVDIS experiments at JLAB as well as electron 
ion collider (EIC) data, may allow a separation of higher twist contributions to Arl from 
the CSV effects depending on their relative sizes. In order to so however, one must have an 
accurate determination of the size of the matrix element of the operator in Eq.(5) as well 
as accurate knowledge of its mixing with other twist-4 operators under the renormalization 
group. 

In light of these experimental developments as well as the broader goal of eludicating 
the transition from perturbative to non-perturbative features of nucleon structure, we have 
undertaken the present study. Our goal in this work is to compute the one-loop anomalous 
dimension matrix for flavor non-singlet twist-4 operators at lowest spin, including all mixing 
effects. Our calculations of the anomalous dimension matrix are an extension of previous 
works [33-35] where parts of the one-loop anomalous dimension matrix were computed. 
In particular, the graphs of Fig. 1(B) (below) were not included in previous analyses, and 
our results for the graphs of Fig. 1(A) differ from those computed in [33]. We defer the 
calculation of the RG evolution of gluonic operators and twist-4 operators of arbitrary spin 
for future work. To this end, we have also listed the iso-singlet flavor operators which mix 
with gluon operators. 

We have combined our perturbative calculations with input from the lattice [36] and 
phenomenological estimates [37] to provide illustrative computations of the evolution of the 
isovector, flavor non-singlet contribution to F2 in a range of Q 2 so as to make contact with 
the CLAS analysis [11]. We demonstrate that theoretically one expects a relatively tiny over- 
all twist-4 contribution to the moment in the resonance region, having a magnitude that is 
consistent with the CLAS analysis. We also show that within twist-4, cancellations or en- 
hancements can occur between different flavor channels contributing to the leading moment 
of F 2 . Thus, a suppression of the twist-4 contribution may be due either to cancellations 
between different operator contributions or to relatively small individual matrix elements 
themselves. Our key results can be summarized in Eq.(56) and Figs. [7] and [8]. These 
results demonstrate that probing higher twist effects in the leading moment of F2 would 
require a substantial improvement in experimental precision. In particular, the observation 
of any breakdown of cancellations due to (Revolution would likely require a substantial 
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reduction in experimental error. As a corollary, we also note that a complete QCD analysis 
of twist-four contributions will require new non-perturbative computations of the twist-4 
operator matrix elements, as the illustrative results given in our study have required making 
an ansatz about the values of several of these matrix elements. 

This paper is organized as follows, in section II we review the standard formalism of the 
operator product expansion (OPE) and establish basic notation. In section III, we review 
the leading twist basis of operators and list the basis of quark and gluonic twist-4 operators. 
In section IV, we present both the Feynman diagrams and renormalization factors for the 
basis of operators introduced in section III, and in section V we introduce a power counting 
scheme to ensure the anomalous dimension has a consistent power in the strong coupling. 
In section VI we organize the basis of twist-4 quark operators in terms of the irreducible 
representations of the flavor group SU(3)^. In sections VIII we list the tree level Wilson 
coefficients used to plot the leading moment of F2 and, in section X, we estimate values for 
the twist-4 reduced matrix elements based on lattice computations and model independent 
estimates. Finally in section IX, we present our results for the leading log evolution of 
F 2 (x,Q 2 ). We discuss these results and comment on future work in section XII. 



II. GENERAL FORMALISM 



In this section we review the formalism and relevant notation for electron-nucleon deep 
inelastic scattering (DIS). The differential cross-section in the one photon exchange approx- 
imation is given by 

d2(7 - a " E L W »v (Q) 

where Q is the laboratory solid angle of the scattered electron, E' is the energy of the 
scattered electron, L^ u is the leptonic tensor 

L^v = 1{k ii k v + k^k u — g^ y k ■ k). (7) 

In L^y, k^ and k'^ denote the initial and final electron momenta respectively and is 
the hadronic tensor given by 

W" u = ImT^, T^ = i J d 4 xe iq - x (N\T(J^x)J u (0))\N), (8) 

where J M is the electromagnetic current of the struck quark and g M = k^ — k' 11 with q 2 = —Q 2 . 
The above form of follows from the optical theorem in which the imaginary part of 
the forward Compton amplitude is related to the cross-section for fully inclusive scattering 
off the initial state nucleon (N). Lorentz and gauge invariance dictate the following general 
form for the hadronic tensor 



(3 



where P M is the initial nucleon momentum, v = P ■ q and are dimensionless structure 
functions. The moments of these structure functions can be written in terms of the OPE as 
shown in Eq.(2) for F 2 . 

The structure of the product of electromagnetic currents given in Eq.(8) at light-like 
distances is given by Wilson's operator product expansion, see Refs. [34, 38, 39] 

00 

J (I) J K) = C U^) G^"""(0) x m (9) 

n=0 i,T 

where we have suppressed the Lorentz indices on the currents. The operators appearing on 
the RHS of Eq.(9) are symmetric and traceless in indices fii . . . fi n and thus have a definite 
twist (dimension - spin) denoted by r. 

The structure functions Fi are determined by both the Wilson coefficients and the matrix 
elements of the operators in Eq.(9), the Fourier transforms of C™ T (x 2 ) are related to the 
moments of the structure functions e.g. 

T dxx n ~ 2 F h2 (x,Q 2 ) ~ ^q r (Q 2 )(iV|Or(0)|iV). (10) 
Jo j 

The dependence on Q 2 is controlled by the anomalous dimension of the operators in 
Eq.(10). The bare (On^) and renormalized (O™) operators of Eq.(2) are related by 

(ryri(b) _ yij /rirj /-. -| \ 

where Z% T denote the renormalization constants. They are, in general, matrices since differ- 
ent operators of a given spin n mix under renormalization. From the scale invariance of the 
bare operators one can derive the RG evolution equations 

A A 

^dn n ~ 7 " 7 ~ nT djA ' W 

and from the //—independence of the moments it follows that the Wilson coefficients satisfy 
the RG equation 



/^ci = T £ c; r) (13) 



which can be solved to give 



cUQ 2 /v 2 ,g) * E c ir(^9(t')) T ex P {- f MinMt')) 



(14) 



JI 



Where t = 1/2 ln(Q 2 //i 2 ), and g(t) is the running coupling in QCD. In what follows, we 
will evaluate 7,^ and its eigenvalues for the non-singlet, twist-4 operators. Our main phe- 
nomenological task will then be an evaluation of Eq.(10) in the resonance region. 
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III. OPERATOR BASIS 



In this section we review the basis of operators that appear at twist- 2 and twist-4. At 
twist-2, it is well known that there are just two towers of operators for a given spin n. The 
multiplicatively renormalizable flavor non-singlet (NS) operators are 



q,fil...fi n 



— trace terms, (15) 
and the two types of flavor-singlet operators that mix under renormalization are 
^Li.-Mn = iN ~ ls $f 7mi^2 • • • D„ n ip f ] - trace terms, 

CS ;w ...Mn = 2 ^^ S [W • • • - trace terms. (16) 

Here ipf denotes a quark field of flavor / and F al3 denotes the gluon field strength tensor and 
A a is an SU(3) / genertor. The operation S reminds one to symmetrize the Lorentz indices in 
brackets. The operators in Eqs.(15) and (16) are thus completely symmetric and traceless in 
the indices ■ ■ ■ \i n and transform under irreducible representations of the Lorentz group of 
spin-n. The anomalous dimension matrix for the flavor singlet operators takes the schematic 
form 

„« - I a ff a 'f9 



7 » = "ti ~t? . (17) 
V a 9 f a 99 J 

The diagonal entries all t and a™ arise from self-renormalization graphs for operators ,, 
and Oq.^...^ respectively. The off-diagonal entries come from graphs that mix these two 
operators. The flavor non-singlet operator 0^ s undergoes multiplicative renormalization 
since it cannot mix into the flavor-singlet operators O^, Oq. 

The situation for twist-4 is more complicated. In general the operators at twist-4 can be 
classified into several types which mix at the one-loop level. The specific number of operators 
grows with the spin-n unlike the case at twist-2. A complete basis of twist-4 operators with 
tree level Wilson coefficients was first given in Refs. [20-22] and parts of the anomalous 
dimension matrix were computed at one loop in Refs. [33-35]. The 'canonical' basis in 
Refs. [20-22] was constructed by requiring that the time-ordered product of electromagnetic 
currents is expanded in terms of operators that (a) are totally symmetric, (b) traceless, and 
(c) contain no contracted derivatives. In this paper, we extend the work of Refs. [33-35] 
and compute the full anomalous dimension matrix of flavor non-singlet operators at twist-4 
and spin-2. The operators at twist-4 that contribute at spin-2 are given by 
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A.gw) = gi, R r a ^rr^ R ^r~ 2 ' k ' l r a ^ 



L ■ 



a . g? ) = ^ r 775*" V, 



(18) 



where A is a light-like vector, A • g n = A^ 1 ■ • • A^Q^...^, d = iA^Du, f = F p ^A p , 
and = e pl3(TT F aT A p . The subscripts i?, L on the quark fields denote the chirality so that 
ipR,L = ^p-ip- I n this paper we compute the anomalous dimension matrix of the operators 
listed in Eq.(18), however in the small x-Bjorken domain, we expect purely gluonic operators 
to make the main contribution. For twist-4, in addition to the quark operators listed above, 
we list purely gluonic operators as well [40]. 





A-0 G1 

n 


= Ti[F a/3 d n 




A- 


n 


— m — 4- 

= Tr[f a d 


-k-i ->fc ->l 

Fd f p d f 


A- 


^ n 


— m — 4- 

= Tr[/ Q d 


/ad 


A- 


f)G4(k,£) 
w n 


— S-n — 4- 

= Tr[/ a d 


/^ //jd n 


A 


• of w 


— >n-2- 

= Wad 


_< i^ d>] 



(19) 

IV. RENORMALIZATION OF TWIST-4 OPERATORS 

The renormalization of the twist-4 operators listed in the last section is complicated by the 
large number of mixings present. In this section, we classify the operator basis into distinct 
types to better organize the calculation of the one-loop anomalous dimension matrix that 
determines the RG evolution at leading order. In Eqs.(18) and (19), the operators Q\ — 
are 4-quark operators, Q 7 n , are 2-quark operators, and Q G1 —Q^ 5 are pure gluon operators 
which we symbolically denote as 4Q, 2Q, and G type operators respectively. In terms of 
this classification, the anomalous dimension matrix then takes the following schematic form 

7 4Q^4Q 7 4Q^2Q ytQ-J-G 

| ll Q ^ AQ ll Q ^ 2Q ll Q ^ G | , 
G^4Q ^G->2Q ^G->G 
In In In 
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(A) (B) 

FIG. 1: Renormalization of 4Q operators. The diagrams in (A) are self-renromalization graphs. 
The graphs in (B) are annihilation graphs. 



where "f^ 4 ^ s a matrix that arises from the self-renormalization graphs of the 4Q opera- 
tors, Jn 4 denotes contributions from graphs that mix the 2Q operators into 4Q operators 
and so on. Recall that we are restricting our analysis to spin-2 (n=2), flavor non-singlet 
twist-4 operators. At one- loop, Fig. 1(A) shows the QCD self-renormalization graphs of the 
4Q type operators. The graphs of Fig. 1(B) contribute to the 4Q— S-4Q self-renormalization 
after using the QCD equations of motion which were not considered in previous work. Fig. 2 
shows the self-renormalization of the 2Q operators, and Fig. 3 shows the 2Q— )4Q mixing 
graphs. For the graphs in Fig. 2, we have chosen to compute using the background field 
method [41]. 

Below we review in schematic notation, the ingredients that go into the anomalous di- 
mension calculation. The bare and renormalized operators are related as 

Q\ = Z ij O j , (21) 

where the subscript b on the LHS indicates a bare operator and the operator on the RHS 
denotes the renormalized operator and Z l i denotes the renormalization constants. The 
indices i,j run over the basis of operators. We also denote the renormalization factors 
for the massless fermion wave function and the strong coupling constant as Z^ and Z g 
respectively so that the bare (b) and renormalized quantities are related as 

ip b = y/Z^ip, g b = lfZ g g. (22) 
The renormalization constants can be expanded around unity as 

Z 4 , = l + 6Z^, Z g = 1 + 6Z g , (23) 



10 



where 5Z^ and SZ g denote the contributions from higher order perturbative diagrams. 

We have found that the anomalous dimension for ^ 4( 3 _i>2< 3 i s zero as one expects on general 
grounds [26]. Consequently, only the 72^"^, 72 < ^~ >2< "S an d j 2 ®^ 4 ® blocks of the anomalous 
dimension matrix are relevant. We break the matrix Z l i of Eq.(21) into the component 
blocks Z 4 ®' 4 ®, Z 2 ®' 2 ®, and Z 2 ®' 4 ® corresponding to mixings among the 4Q operators, the 
2Q operators, and the mixing of 2Q operators into 4Q operators respectively. As mentioned 
previously, the 4Q operators do not mix into the 2Q operators and since we restrict our 
analysis to flavor non-singlet operators we do not include the pure gluon G-type operators 
in the basis. The Z 4 ®' 4 ® and Z 2 ®' 2 ® renormalization matrices can be expanded around the 
unit matrix as 

Z 2Q,2Q = 1 + SZ 2Q,2 Qj Z 4QAQ = j + 5Z *QAQ ^ ( 24 ) 

while the off-diagonal block Z 2 ®' 4 ® gets non-zero contributions starting at one-loop and is 
written as 

Z 2QAQ = 5Z 2Q,4Q ( 25) 

We now have all the necessary notation to discuss the extraction of the one-loop anomalous 
dimension. We outline the steps for the 4Q and 2Q operator renormalization below. 



A. Four-Quark Operators 

The bare 4Q operators have the schematic form 

where we have suppressed flavor indices and the Lorentz and Dirac structure. The renor- 
malized and bare 4Q operators are related as 

Q 4Q = [Z~ 1 ) 4Q ' 4Q o 4Q 

= Tpi/j^-ip + (25 Zj, + 25Z g - 5Z 4QAQ ) # V e Tptjjipip, 

(26) 

where the second term above is just the counterterm and determines the renormalization 
matrix 5Z 4 ®' 4 ®. Once 5Z 4 ®' 4 ® is extracted from the counterterm above, the anomalous 
dimension matrix is given by is given by 

ytQ^4Q _ ^QAQyl^^_ z iQAQ (27) 
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FIG. 2: Feynman diagrams for the renormalization of 2Q operators. 



B. 2Q Operators 



The bare 2Q operators have the schematic form 



Ol Q = i> b g h F h if) h 



(28) 



where Ff, denotes the bare field strength tensor and we have suppressed flavor indices and 
Lorentz and Dirac structure. There renormalized 2Q operator is related to the bare operators 
as 

= ijgF^ + (5Z^ - 5Z 2Q > 2Q ) iigFij + (5Z^) 2Q ' 4Q Z^^g 2 ^T^T^, (29) 

where the two terms in the first line above correspond to mixing among the 2Q operators 
and the mixing of the 2Q operators into 4Q operators respectively. The combination g^Fb 
remains unrenormalized in the background field method, and the last two terms in the 
second line of Eq.(29) denote the counterterms and the anomalous dimension components 
are given by 



7 



dfi 



(30) 



For one-loop renormalization, Eq.(29) simplifies to 



2Q = ^gF^j + {8zf - 6Z^ 2Q > 2Q )4,gF^j + (^m^^j-WW^, 

where the superscript (1) on the renormalization constants indicate the respective one- loop 
contributions. 



12 



FIG. 3: Feynman diagrams for 2Q — > AQ mixing. 
V. CONSISTENT POWER COUNTING IN POWERS OF g 

Before proceeding to the calculation of the anomalous dimension, we address an issue 
concerning consistent treatment of orders in perturbation theory [20]. Using abbreviated 
notation, we collectively call Oi = {ipip) 2 , 2 = ipFip, G 3 = FD 2 F, and G 4 = F 3 . A close 
look at the mixings among 0\ — G4 reveals that Zg 3 ^g 4 is order g 3 whereas Zg^g z is order 
g as shown in Fig. 4. One can readily see that the counting inconsistencies persist when 
computing the mixings Zg s ->Oi as well. It is desirable to write the anomalous dimension in a 
scaled form 7^ ~ g 2 <% when computing Eq.(10). This form of 7 renders direct calculations 
of the integral 



T exp 



W) 



(31) 



to be quite simple. However, a leading log evolution of the operators 0\ — G4 leads to 
an anomalous dimension matrix which is not proportional to one consistent power in the 
coupling, e.g. 



I 







\ 



9 2r Tu 9 3r Ti2 

9121 g 2 122 gi23 g 2 12A 
g 3 i32 g 2 i33 g 3 i3A 
g 2 742 9 743 g 2 744 / 







(32) 



A form of the mixing matrix proportional to g 2 in lowest order can be regained by an 
appropriate rescaling of the twist-4 operators. We have chosen the following redefinitions 



2 g{i>^), 
G 3 -> FD 2 F, 
G 4 -> gF 3 . 
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FIG. 4: Left - A representative Feynman diagram illustrating operator G4 mixing into G3, an 
order g correction to G3. Right - The Feynman diagram illustrating the mixing of G3 into G4, 
which is order g 3 . 



Of course, the dominant logarithm is independent of such conventions. After these redefini- 
tions, the anomalous dimension matrix has a homogenous scaling in the strong coupling 



7(s') ~ £ 



( 7n 7i2 \ 

721 722 723 724 

732 733 734 

\ 742 743 744 / 



(33) 



VI. FLAVOR STRUCTURE 



In this section we discuss and establish notation for the flavor structure of the twist-4 
operators. The structure of the anomalous dimension matrix can be organized according to 
flavor structure since QCD with massless quarks preserves flavor symmetry. The electro- 
magnetic current entering in the forward Compton amplitude in Eq.(8) is given by 

J"(x) = ^(x)YQ^(x), Q = i(A 3 + i=A 8 ), (34) 

where ip is a column vector in flavor space so that ip = (ip u ,ipd,ips) and Q is the electro- 
magnetic charge operator acting on ip and can be written in terms of the SU(3) / Gell Mann 
matrices Aj as shown. The twist-4 operators from the OPE of the product of electromagnetic 
currents in Eq.(8) can be classified in terms of their transformation properties under SU(3) /. 
Schematically, the 4Q and 2Q operators have the following flavor structures 

4Q : A) rjjQWQij), B) ipQ 2 ipipip, 

2Q : C) ^Q 2 ^, (35) 

where the precise color and Dirac structure is suppressed. The flavor structure A) in Eq.(35) 
arises from the first handbag diagram of Fig. 5. The second diagram of Fig. 5 generates 
both B) and C) flavor structures where the flavor structure in B) arises after an application 
the gluon equation of motion (EOM) for the external gluon. These flavor structures can 
then be decomposed into irreducible representations of SU(3)y with definite isospin (I,I Z ) 
and hypercharge (Y=2Ag/v / 3)[42]. Since the charge operator Q is a linear combination of 
A3 and Ag all these operators have I z = Y = 0. 
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FIG. 5: Left - Double hand bad diagram with flavor structure ipQipipQip. Middle - Feynman 
diagram with flavor structure ipQ 2 ip. Right - Feynman diagram with flavor structure ipQ 2 ijjil>xjj 

The flavor decomposition of the 4Q operator of type A) in Eq.(35) is given by 

#w WW = \P~ o 2 /4 + ^=o 2 /4 + -%=of4 + -?=o& + -^=oKi - — o}4 

y 3 /-a ^ /-i /-o /-i 7-0 3 7-0. 



(36) 

and the flavor structure ijjQ 2 ip that appears in the 4Q operators of type B) and 2Q operators 
of type C) in Eq.(35) is decomposed as 



4>Q 2 ^ = 0%* + O?f 



(37) 



The superscripts on the operators on the RHS of Eqs.(36) and (37) denote the SU(3)/ 
representation and the subscripts denote the isospin representation of the SU(2)j subalgebra 
of SU(3) /. The labels (A, B, C) are included to remind the reader the specific flavor structure 
given by Eq.(35). For notational convenience we define "meson" fields that make the flavor 
structure of quark bilinears manifest as 



7T 



7T 



7T 



(ill: 
l 



K + = su K° = ds Vs = -j=(dd + uu — 2ss 



V2 
—ud 



(dd — uu) 



K° = sd 



K~ 



-US T] 



'uu + dd + ss) 



In terms of these fields, the flavor structure of the operators of the type A operators appearing 
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on the RHS of Eq.(36) is given by 
v6 

OT = f ~ -^=[{K°K° + K~K + ) + (K°K° + K + K~) + ^ As + m*°)] 
v 10 



Of J ~ 4= [ g (*V - TT+vr- - 7T-7T+) - - K-K+) - \{K»K» - K + K~) + 



Of* ~ - ^1 [ ^i, + + + K + K~) - -L(As + ^°)] 

O^ ^[(vr°7T - 7T+7T- - 7T7T+) - -(K°K° - K~ K + ) - -(K°K° - K + K~) - W8 ] 

y5 2 2 

O)f U(?rV - 7T + 7T- - 7T-7T + ) + - + ) + (K° K° ~ K + + ^g] . 



(3* 



The last operator in Eq.(38) is a flavor singlet and can mix with the pure gluon 

operators and will not be considered in the rest of the analysis. The flavor structure of 
operators of type B in Eq.(37) are given by 

°f=o -^=W+3^V, (39) 

and the flavor structure of the two quark operator ipQ 2 ip is, 

05 ~ ^ " 171* (40) 



VII. STRUCTURE OF THE ANOMALOUS DIMENSION MATRIX 



In this section, we expand the discussion of section IV on the structure of the anomalous 
dimension matrix to incorporate the flavor structure discussed in section VI. Equations (36- 
40) give the SU(3) flavor decomposition of the 4Q and 2Q type operators. The conservation 
of flavor in massless QCD implies that the Of^i o operators in Eq.(36) will not mix with 
operators living in a different representation of SU(3) or with those in a different isospin 
subgroup. On the other hand, the octet operators I=1 and Oj'_{ can mix with each other. 
Thus, the analog of Eq.(21) that relates the bare and renormalized operators for OjL 210 
(dropping the A label) is diagonal 

of =2 \ /P /=2 
6f =l = o p J=1 o 
6f = J h \ o o p /= 









oil, 







(41) 
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where the vector Of is a 6- dimensional column vector corresponding to the Dirac and color 
structures of the 4Q operators QfcZ-f °f Eq.(18) with flavor structure given by the 27 flavor 
representation with isospin I appearing in Eq.(36) 

n 2 ? - (n 1 ^ n 2 (°'°) n 3 (°'°) n 4 (°>°) n 5 (°<°) n 6( -°>°h T (ao\ 

W I — lVn=2 > Vn=2 iVn=2 iVn=2 ) Vn=2 iVn=2 J27,/' l^ Z J 

where the superscript T denotes the transpose. Note that there are no 2Q operators in the 
27 representation of SU(3) /. The renormalization constants P/ are thus 6x6 matrices. The 
anomalous dimension matrix for the 27 operators is given by 

■ff = P 7 V|;P, (43) 

The octet sector of the anomalous dimension matrix is more complicated. Operator types 
(A, B,C) in Eq.(35) all contain flavor-octet operators, and in general, these three operator 
types will mix under renormalization. For flavor structures of type A) and B) one encounters 
only four quark operators while for C) one has the two-quark operators. For convenience, 
we embed those of type A) in a ten-component vector Of A and combine those of type B) 
and C) into a second ten component vector, 8,BC . The first 6 entries of 8 ' BC are filled by 
the Q 8,B (4Q) operators and the last four operators are filled by Q 8,c (2Q) operators, 

A8,A _ / n l(0,0) n 2(0,0) n 3(0,0) n 4(0,0) n 5(0,0) ^6(0,0) n n n n\T 
W I — Wn=2 > Vn=2 > Vn=2 ) Vn=2 ) Vn=2 ) Vn=2 ) u ) u ) u ) u )8,A,Ii 

n^BC _ / n l(0,0) n 2(0,0) n 3(0,0) n 4(0,0) ^5(0,0) ^6(0,0) ^7(0) ^7(1) ^8(0) n 8(l) xT 

W I — Wn=2 i V n =2 i V n =2 > V n =2 i V n =2 > Vn=2 > V n =2' V n =2' V n =2' ^n=2J8,BC,I- 

(44) 



With these definitions, the bare and renormalized operators are related as 
The matrices L, M, N, Q are then 10 x 10 matrices which have the form 

jAQ^iQ jAQ^2Q \ / |^4Q^4Q |^4Q^2Q 




(45) 



j^2Q->4Q j^2Q^2Q J ^J=l — I ^2Q->4Q ^2Q^2Q 

y^Q^AQ y^Q^2Q \ / q4Q^4Q q4Q^2Q 

M 2Q^4Q M 2Q^2Q J Ql = ( q2Q->4Q q2Q->2Q 



(46) 



The Lf and Nj matrices encode the remormalization structure of the 4Q and 2Q operator 
structures for the 8AB and 8C representations respectively. The matrices Mj and Qj encode 
the mixing of the 4Q and 2Q operator structures between the 8AB and 8C representations. 
Many of the submatrices in Lj, M/, N/, Qi vanish 

^4Q^2Q = jJQ^Q _ fJQ^Q — M 4 j Q ^ 2Q = Nj Q ^ 2Q = Q 4Q ^ 2Q = (47) 

All but the h 4 j Q ^ 2Q vanish since only AQ operators appear in 8 j' A . The remaining sub-blocks 
M^ - * 2 ^, p^Q^ 2 ^ q4Q->2Q wou j^ gj ve r j ge £ mixing of AQ into 2Q operators (see Section A 
for more details). However, it is known on general grounds that such mixing does not arise 
(see e.g., Ref. [26]), a result that we have reproduced with our explicit computation. 
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VIII. WILSON COEFFICIENTS FOR n = 2 



Having outlined the flavor structures of the twist-4 contributions, we now present the 
Wilson Coefficients for the twist-4 2Q and 4Q operators at leading spin. According to the 
formalism established in Ref. [20], the Compton amplitude at twist-4 naturally divides into 
two pieces arising from the graphs of Fig. 5 



dW^T[J M (x)J,(0)] = X^ + Yi 



(48) 



The Y^v term arises from the double handbag-type diagrams (corresponding to the first and 
last diagram in Fig. 5) and the X^ v term arises from the remaining diagrams. The explicit 
calculations of X^ v and Y^ are given in detail in Ref. [20], and we summarize the full form of 
these expressions in appendix B. For a leading moment (n = 2) analysis, these expressions 
simplify greatly. The Y^ v term is given by 



T=4,n=2 

/IV 



(IV 



M r^ 2 C9 2(0 ' 0) 

1 fj,v °n=2,^i/i2 



1-12 



(49) 



and the Xn U term is given by 



X 



T=4,n=2 



g_ 

g 

2g 6 



fJ,U 



IL°»-3 ? .0£>-3 ? .0S} 



p-q 



g 



(p-q) 2 . 



\q-OTi + \q-0^ + 5 q .O^}. 

(50) 



Here q ■ O is shorthand for g Ml . . . q^O^ 1 "' 11 " , while the explicit form of the operators ap- 
pearing in X^ y and Y^ u in terms of the canonical operators of Eq.(18) is [20] 



A- 


w n=2 


= A-QlT - 2 A-QS 0) 


+ a-qIT 


A 


n m 

' w n=2 






A 


' w n=2 


- A-QS 




A 


' w n=2 


= A-Q« + 2A-g^ 





(51) 

Note that the first term in the RHS of Eq. (50) contributes to the longitudinal structure 
function Fl while all other terms in Eqs. (49,50) contribute to F 2 whose moments we analyze 
in this work. 

The foregoing decomposition does not yet reflect any flavor structure. To avoid introduc- 
ing unwieldy notation, we simply indicate this structure below: 



w n=2 



A-O 



7(0) 



A-O 



n=2 

3(0,1) 
n=2 



Flavor Structure A: 27,8 
Flavor Structure B: 8 
Flavor Structure C: 8 . 



(52) 
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We now apply these results to the forward matrix elements of the vector current correlator 
in Eq.(48). Several steps are required: (i) expressing the matrix elements in terms of those of 
the individual twist-four operators Q % ^2\ (ii) decomposing the current-current product in 
terms of the various flavor structures A), B), and C); and (iii) expressing the latter in terms 
of the operators associated with their SU(3) flavor decomposition given in Eqs. (36,37). 
Starting with the first of these steps, we write the matrix elements of the twist-4 canonical 
basis operators as 

(N| Q£W> IN) = A*™ ( Pltl p„ - \M 2 N g^ (53) 

where the A factors are reduced matrix elements encoding the non-perturbative multi-parton 
correlations. Second, using this form for the matrix elements in Eqs. (49) and (50), we write 
the n = 2 component of T MJ , as 

= -i ! d 4 xe^(P|T [^(2)^(0)] \P) 

J ™= 2 

w 2 d Ui/ f . 5 . 1,1 cu 2 e uu ( 1 . 3,1 

= —V- \ Aa + s Ae + lei M " U A ° " 8-M ■ (54) 

in agreement with Ref. [37]. Here, we have defined u 2 = l/x 2 B and have indicated the flavor 
structures of each matrix element for clarity, introducing the shorthand notation: 

A A = A Um - 2 ^ 4 ' (0 ' 0) + ^ 6 ' (0 ' 0) Flavor Structure: $Q^Q*/j 

A B = A 2 ' m + 2 ^ 4 ' (0 ' 0) + ^ 6 ' (0 ' 0) Flavor Structure: 4>Q 2 ^ 

A c = A 7 ' {0) - A 7 ' {1) Flavor Structure: ipQ 2 ifj. (55) 

The tensors e^ v and d^ v are written in full in section B. The coefficients of d^ u contribute only 
to the F 2 structure function whereas the coefficients of e^ u contribute to the F L structure 
function [20]. 

We now express the leading moment of the isovector part of F 2 in terms of the foregoing 
matrix elements. In doing so, we also carry out the SU(3) decomposition following [33]. The 
result is 

Ml2l T=4 (Q 2 ) = Jdx B F^U^Q 2 )= ^E{^ C A j (Q 2 )^l 

+ -^=0^ (g 2 )^- + C 8 B \Q 2 ) A%, + C S C \Q 2 ) A c ,) (56) 

= M 27 (Q 2 ) + M 8A (Q 2 ) + M 8B (Q 2 ) + M SC (Q 2 ). (57) 

We have introduced the notation A for the matrix element of Of =x appearing in Eq.(38). 
The subscripts A, B, C indicate the specific flavor structure of each operator, e.g. that Ab 
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TABLE I: Tree-level Wilson Coefficients C^^(Qq) evaluated at the input scale Qq. Here, ./V denotes 
the SU(3) multiplet while j runs over the set of isovector canonical operators in Eq.(18). We have 
not included Wilson Coefficients for operators with j = 8(0) and j = 8(1) as these coefficients all 
vanish at the input scale.) 





i = l 


3=2 


J = 3 


i = 4 


j = 5 


j = 6 


3 = 7(0) 


3 = 7(1) 







1 





-2 





1 













1 





-2 





1 













5/8 





5/4 





5/8 




























1/16 


-1/16 



is the matrix element of an operator of type ipQ 2 ^^^ and Ac is the matrix element of type 
^jjQ 2 ^p respectively 1 . The explicit factors of 2/y/lO result from the SU(3) decomposition of 
the A, B, C type operators while the C^ 7j '(<5 2 ) are the corresponding Wilson coefficients. 
The index j runs over all relevant operators in the canonical basis in Eq.(18) [20]. The 
values of the Wilson coefficients at an appropriate input scale (discussed below) are given 
in Table I. 

One may ask whether it is possible to isolate experiementally the 27-plet and octet con- 
tributions to the leading moment. To this end, we note that for unpolarized D1S processes, 
the photon couples to a current we denote by J7~em — Vs + where V signifies the 

vector nature of the current, and we have defined Vf" = ■0Tpy'V- Following [33], a similar 
flavor decomposition can be done for the electroweak charged current. The charged current 
contains both a strangeness-changing piece and a non-strangeness-changing piece. For the 
strangeness-changing part of the charged current, the isovector 27-plet and octet moments 
are expressed as linear combinations of moments of F 2 extracted from neutral and charged 
current DIS processes (see Ref. [33] for details) so that we can write 

M*l 1>n=2 = JI. h{M? ~ M?) - ^V(M 2 yp + Mf - M 2 - - Mf ) ) , (58) 
V 10 ^ 2 sin Vc J 

and 

MjU, n=2 = ^5 {2(M 2 cp - M?) + -ij-W + Mp - M 2 - - Mf )) . (59) 
10 [ 2 sm Vc ) 

Thus, through a combination of experiments it is in principle possible to isolate specific 
flavor structures that contribute at twist-4 to the moments of F 2 . 



1 For the flavor 27, the subscript A is clearly redundant, as only the structure A can yield a 27 after SU(3) / 
decomposition. However, we retain the subscript in this case for overall uniformity of notation. 
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IX. LEADING LOG RG EVOLUTION 



In this section we present results for the Wilson coefficients of the flavor non-singlet twist- 
4 operator combinations. Within each flavor representation, the RG evolution is affected 
by mixing between various operators in the canonical basis. The evolution of the Wilson 
coefficients in flavor representation R is then given by 



/' 2 



■9(t) 



3 



Cf(l,(?(0))Exp 



dt'lW)] 



j i 



(60) 



where t = l/21n(Q 2 //U 2 ) and the subscripts i,j label the Wilson coefficients of the canonical 
operators in Eq.(18). The Wilson coefficients (7^(1,^(0)) correspond to the values obtained 
in the matching calculation at /x 2 = Q 2 . For leading log running, the C^(l, g(0)) correspond 
to the tree level values obtained from the OPE. The Wilson coefficient on the LHS of Eq.(60) 
corresponds to the value of the Wilson coefficient after RG evolution from the initial scale 
Q 2 to the final scale /i 2 . With one loop running the evolution of the strong coupling is given 
by 



<? 2 (0) 



2/W(0)t' 



(61) 



where (3 = l/(47r) 2 (ll/3 Ca — 4/37/ n/) and g(0) corresponds to the strong coupling evalu- 
ated at fi 2 = Q 2 . To solve the evolution equation, we first diagonalize the one- loop anomalous 
dimension matrix so that 



lji[9{t)} = 9 2 {t) Rjm d m l R t >, 



(62) 



where d m e = S m £d m is the diagonalized matrix with eigenvalues d m and R denotes the 
appropriate rotation matrix. The anomalous dimension matrix elements 7^ for the various 
canonical operators in different flavor representations are presented in section VII. The 
evolution equation for the Wilson coefficients can now be written as 



a 



R 



^,g(t) 



E Cf(l,p(0))Exp 



3,m 



g 2 (0)d r , 



dt' 

l + 2/W(0)t' 



The Wilson coefficients are evolved from the scale Q 2 to the scale of the non-perturbative 
matrix elements which we denote as /1 2 = Ql and refer to as the input scale. Thus, in terms 
of the latter the Q 2 dependence of the Wilson coefficients is given by 



E ef(i,<?(o))Ex P 



3,m 



t0 dt ' g2[ V )dm 

l + 2/3 ^ 2 (0)t'J 



RjmR m ], (64) 



where to = l/21n(Q /Qq) and (1,(7(0)) are the tree-level values determined from the 
OPE at the matching scale Ql = Q 2 . 
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FIG. 6: Evolution of Wilson coefficients using Eq.(67). The flavor representation of each coefficient 
is indicated at the top of each plot (see Eq.(56)). As input to Eq.(67), we have taken A = 0.3 GeV 
and rif = 3. 



The above expression can be simplified further to give 
'Q 2 



where 



so that 



Ql 



2 >9{to] 



J2 C?(l,g(0)) l + /3 /(0)ln 



3,m 



QV 



dm 
"2/3 



c 



R 



Ql 



X)cf(i^(«)) 



3,m 



HQVQl) 

ln(g 2 /A 2 ) 



dm 
"2/3 



D D — 1 



(65) 



(66) 



(67) 



In Fig. 6, we plot the RG evolution of the Wilson coefficients in the isovector combination. 
We have used the input scale Ql = 5 GeV 2 , corresponding to the scale at which existing 
theoretical evaluations of the non-perturbative matrix elements have been performed (see 
Section X below) and the tree-level Wilson coefficients given in Table I. We have chosen 
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rif = 3 in Fig. 6 which is consistent with the SU(3)/ decomposition. At the top of each 
plot, the notation C^ BC denotes the Wilson coefficient in representation R and of the type 
A,B,C corresponding to flavor structures ipQipipQip, ipQ'^ipipif), and ipQ 2 ip respectively. 
Note that we evolve the coefficients to smaller rather than larger values of Q 2 since we are 
interested in the kinematic regime relevant to the CLAS analysis and since the scale Qq 
associated with existing theoretical matrix element input lies above this region. 

Several features emerge from Fig. 6. First, several of the Wilson coefficients that vanish 
at the input scale become non- vanishing and significant in magnitude at lower scales due to 
operator mixing. Second, the coefficient C^f for j = 7(1) changes sign, going through zero 
near Q 2 = 2 GeV 2 . All other coefficients increase in magnitude with decreasing Q 2 . We find 
no satisfying explanation for the particular behavior of for j = 7(1) other than that 
its tree-level value is rather small and that there exist various contributions from operator 
mixing having opposite signs. Together, these features lead to the possibility that significant 
cancellations among various operator contributions may occur at one scale - leading to a 
suppressed twist-four effect - but that these cancellations are broken at other scales by the 
Q 2 evolution. Alternately, the overall twist-four contribution to the leading moment may 
be relatively small at all scales due to suppressed values of the operator matrix elements at 
the input scale. We explore these possibilities in the following sections. 

X. HIGHER TWIST MATRIX ELEMENTS - THEORETICAL INPUT 

In order to obtain a prediction for M^zl T=A (Q 2 ), we now require values of the hadronic 
matrix elements of the twist-four operators at the input scale Q (see Eq.(56)). Knowledge 
of higher twist quark and gluon correlators in hadrons is of fundamental interest in order 
to understand the structure of baryons and mesons on the basis of QCD. These matrix 
elements are between hadronic states, making model independent computations challenging. 
However, there do exist attempts in the literature to compute twist-4 matrix elements at 
leading spin on the lattice [36]. Due to the complicated mixings with lower dimensional 
operators however, the analysis performed in [36] was restricted to operators of the specific 
flavor channels outlined in section VI. An alternate phenomenological approach was used 
in Ref. [37]. In what follows, we draw on the results from these two studies to determine 
non-perturbative input for an initial, illustrative analysis of M^zl T=4 (Q 2 )- 

Before proceeding we observe that the aforementioned matrix element calculations pro- 
vide only partial input for the evaluation of the RHS of Eq.(56). The lattice study of Ref. 
[36] gives only a value for the combination A 2 ^ appearing in Eq.(55), at the input scale, while 
the work of Ref. [37] gives only the octet contributions A 8 B c . Away from the input scale, 
one no longer has the specific linear combinations given in Eq.(55), owing to the evolution 
of the Wilson coefficients. A robust prediction for the evolution of M^z^ r=4 (Q 2 ) would 
require values for the individual contributions to Aa,b,c- Moreover, we could find no eval- 



23 



uation of the octet contribution A\, nor do there appear to exist any computations of the 
matrix elements whose Wilson coefficients vanish at the input scale but become non-zero at 
lower scales. Consequently, we will adopt some reasonable ansatz for the individual matrix 
elements, motivated by the computations that do exist but with the caveat that a complete 
computation will require values for all matrix elements. 

To proceed, we first consider the combination A 2 J. The 4Q operators introduced in 
Ref. [36] have been evaluated on the lattice using Wilson fermions in the quenched approxi- 
mation. Phenomenological constraints for twist-4 2Q operators were estimated in Ref. [37]. 
The reduced matrix element of the 4Q operator A 2 J shown in Eq.(56) at an input scale of 
Qq ~ 5 GeV 2 was computed on the lattice and found to be 



A 



27,1=1 



Latt 



■10.4 ± 1.6) x 10^ 4 GeV 2 . (68) 



We take this to be the value of the reduced matrix element A 2 J appearing in Eq.(56). 
Neglecting the Q 2 -dependence associated with evolution, the corresponding contribution to 
the Nachtmann moment for F 2 at leading spin is then [36] 



M 2 n 7 ^U (Q 2 ) = £ dx F%£\x, Q 2 ) = -0.001(1)^. (69 ) 

We note that the matrix element of A 2 ^ l=1 causes the moment to be quite small relative 
to the flavor non-singlet twist-2 operator, whose corresponding contribution is 0.14 at the 
same input scale (see Refs.[43, 44]). 

We could find no corresponding lattice computation of A\ which appears in Eq.(56). As a 
benchmark, we first take its value to be equal to A 2 J at the input scale Ql ~ 5 GeV 2 , though a 
larger magnitude for A\ is possible (see below). This leaves the twist-4, 4Q operator of type 
B) and the 2Q operator of type C). We were also unable to find lattice calculations of these 
operator matrix elements in the literature. As an alternative, we use the phenomenological 
estimates obtained in Ref. [37] that rely, in part, on information extracted from DIS data 
from both CERN and SLAG The 2Q and 4Q operators studied in Ref. [37] have the forms 

Q% = ig i> [D^ F ua } 7 Q 7 5 QV, (70) 

Q% = g 2 4>i,Q 2 r a ij^T^. (71) 

After integrating by parts, the 2Q operator Q 9 is the tree-level 2Q combination for the octet 
combination appearing in Eq.(54): 

Q% = ig${p it F va + F va D^'y a 1 sQ^ 



-> ig 



QT-QT- (72) 
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Following [37], we denote the reduced matrix elements these operators A 9 and A 2 , respec- 
tively. For the proton (neutron) they are written, 



= QlK( d) + Q 2 d K 9 diu) , (73) 

where 

K u ( d) = ^(P\u{D + ,F +tM }n 5 u\P), (75) 

K l = Jp( p \(ui + r a u) (u 1+ r a u)\P), (76) 

K l = ^ p \(dl + r a d) (d 1+ r a d)\P), (77) 

and 7 + = l/\/2(7o + 73)- The neutron matrix elements can be obtained from isospin 
symmetry and the strangeness contribution has been neglected to simplify the analysis. The 
range of possible values for K 9 and K 2 (see Ref. [37] for details) are 

- 0.585 GeV 2 < K 9 < -0.238 GeV 2 (78) 
-0.318 GeV 2 < Kl < 0.203 GeV 2 . (79) 

The values for K^' were then computed by introducing a flavor ansatz, e.g. 

K g,2 /Rg ,2 „ I dx(d(x)+d(x))x 

Where u(x), d(x) are the twist-2 parton distribution functions and (3 = 0.476 at Q 2 = 5GeV 2 . 
Finally, the isovector combination corresponds to the difference A 9 ; 2 — A 9 ^ 2 the final results 
are, 

A% = Al-A 2 n = \{l-f})Kl 

A% = A 9 p -A 9 n = l -(l-P)Kl (81) 

from the range of K 9,2 , the final values for the reduced matrix elements appearing in Eq.(56) 
are 

- 0.10 GeV 2 < A% < -0.04 GeV 2 

-0.06 GeV 2 < A% < 0.04 GeV 2 . (82) 

For Aq, we take the central value in this range. However for ^4^, it is possible to choose a 
range of values both positive and negative and have chosen three representative cases, Ab = 
—0.06 GeV 2 , Ab = 0.0 GeV 2 and Ab = 0.04 GeV 2 . We assume a theoretical uncertainty 
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TABLE II: Reduced matrix elements used in computing M^ = \ r= 4 (Q 2 ) • Note that the input for 
the A\ represent two different ansatz. 



Operator 


Flavor Rep. 


Value 


Method 


A 2 J 




-10.4 x 10~ 4 GeV 2 


Lattice [36] 


A 


tpQi/jtpQi/j 


-10.4 x 10" 4 GeV 2 


A 8 A ~ A A 7 \ Lat t 


A 


ipQiftipQij} 


-10.4 x 10~ 3 GeV 2 


A\ ~ 10 x A A 7 \ Lat t 


A 




-0.06 < A B < 0.04 GeV 2 


Phenomenology [37] 


A 


ipQ 2 ^ 


-0.1 < A c < -0.04 GeV 2 


Phenomenology [37] 



in each value to be of the order of the lattice uncertainty for the 4Q matrix element A 2 A in 
Eq.(68). A summary of all matrix element inputs is given in Table II. 

We now observe that in order to predict the logarithmic corrections to the 1/Q 2 scaling 
of the M^z\ T=i {Q 2 )-, we require knowledge of individual operator matrix elements that 
contribute to A 2 A 7,S and ^4^ c since the Wilson coefficients for the contributing operators 
(labelled "j" ) do not have identical evolution. The matrix elements of these operators at the 
input scale (see Eq.(55)) are constrained by the values listed in Table II. These constraints 
are insufficient to compute the individual values of the matrix elements appearing on the 
RHS of Eq.(55) however. To do so and for purposes of illustration, we have assumed several 
relationships among the A j ^ k ' l \ To illustrate our method, we introduce a shorthand for the 
four-quark operators of type A and B appearing in Eq.(55) 

A 3 £ = A 3 ' 2 ± 2A rA + A j > 6 (83) 
A = A 4 ' 1 ± 2A r ' 3 + A 5 , (84) 

where j denotes the flavor structures of the 4Q operators and takes the value A or B (see 
Eq.(55)). The integers (1 — 6) specify the individual matrix element of the basis operators 
appearing in Eq.(18), and c denotes the presence of an SU(3) color generator. Using the 
above relationships, and the computation of the A± c matrix elements from the lattice, we 
have the following constraints for the 4Q operators, 

A A > 4 = ^(Ai' c -A^' c ) (85) 

A Afi = ^(A Ac + A Ac ) - A A ' 2 . (86) 

Thus, A AA is determined by the lattice values for A^ c and A^ ,c , while A A ' 2 remains uncon- 
strained. We take the value of A A ' 2 to be 25% smaller than the matrix elements A± c and 
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use a similar set of relations for the A± matrix elements. The above procedure is applied 
to all 4Q matrix elements regardless of flavor. The values of the octet matrix elements of 
type A (A 8 appearing in Table II) are taken to equal the flavor 27 matrix elements. 

For the two-quark operators, the analysis in [37] along with Eq.(72) constrains the com- 
bination, 

A c = A i,k=v _ jj,k=\ (87 ) 
We have made the following ansatz to compute the individual matrix elements (A 7 ' k=0,1 ), 

A 7,k=0 = _ A 7,k=l = l^C (g8) 

This leaves the 2Q matrix elements A 8 ' k=0,1 appearing in Eq.(18). We have assumed these 
matrix elements are equal to A 7,k=0,1 . A complete Table of each matrix element, their flavor 
structure and the methods used to compute them can be found in Section E. 

XI. LEADING MOMENT ANALYSIS 

Using the input discussed above, we show illustrative results for the 27 and 8 contribu- 
tions for values of Q 2 in the resonance region in Fig. 7. Under the assumptions used for 
the matrix elements, the twist-4 contributions to the moment are typically quite small in 
the resonance region. For A 8 A ~ -A^Latt, the flavor octet contribution to the moment is 
determined largely by the Ab,c matrix elements since .A^Latt is quite small. Moreover, 
for positive values of A%, cancellations occur within the octet sector, further reducing the 
overall octet contribution. 

The total contribution to M^=\ T=i {Q 2 ) is given in Fig. 8. We plot two representative 
cases for differing values of A\. In the left-hand plot, we have taken A 8 A ~ A/^lLatt- In 
the event that the behaviors of the leading proton and isovector moments are similar, this 
value would be disfavored by the CLAS data. CLAS has measured the twist four moment 
to be positive (Ref [11]) in the range of Q 2 shown in Fig. 8. In the right-hand plot, we have 
taken A\ ~ 10 x A/4 |Latt- Here, a range of small positive values of A% yields an isovector 
moment which is qualitatively similar to the CLAS data for the proton. It is possible, then, 
that the CLAS data provide a hint of a hierarchy among the twist-four operator matrix 
elements, though a definitive statement will require data on F 2 for the neutron as well as 
an improvement in overall experimental precision. At present, the experimental error on 
the determination of the twist-four contribution to the leading moment of the proton F 2 
structure function is ±0.015 for Q 2 m 1 GeV 2 [4, 11], while the magnitude of the leading 
twist contribution to the isovector moment is ^0.1 throughout the indicated region. The 
systematic error in obtaining the moments from data consists of genuine uncertainties in the 
data, as well as uncertainties in the evaluation procedure. 
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FIG. 7: Flavor decomposition of M^l\ T=i (Q 2 ), using the logarithmic evolution discussed in Sec- 
tions VIII and IX and the matrix elements in Table IV. The values of the Ab octet matrix elements 
are shown above each plot. The notation (Q%) follows the notation in Table IV. The specific values 
listed above each plot are used for all type-B matrix elements contributing to the moment. Central 
values were used for the two-quark octet operators for each curve. 



For comparison to the theory predictions, we note that the CLAS data for F 2 is extracted 
from electron-proton DIS and thus is not the linear combination of moments shown in 
Eqs. (58,59). The experimental data are also not separated into flavor singlet and non-singlet 
channels, and thus include the effects of the gluonic operators shown in Eq.(19). One can 
see from Fig. 8 that the experimental precision in extracting the leading moment for higher 
twist does not yet allow one to disentangle the separate flavor channels computed here. Next 
generation experiments may be required to probe the higher twist flavor contributions to 
the leading moment of F 2 . 

Finally, we explore the possibility that the twist-4 contribution may be suppressed at one 
scale due to cancellations between operator contributions and that Q 2 -evolution may lead to 
a breakdown of this cancellation at other scales. In Fig. 9, we have tuned the values of A% 
and Aq such that the total octet contribution to the isovector moment cancels the flavor 27 
contribution at Q 2 w 2 GeV 2 . The total moment is given by the dashed curve in Fig. 9. The 
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FIG. 8: Total contribution to M^ T=4 (Q 2 ), see Eq.(56). In the left-hand plot we have taken 
A\ ~ A 27 \ Latt> i n the right-hand plot we have taken A\ ~ 10 x -A^Latt- I n each plot, the three 
curves represent different values for the reduced matrix element A%. We have taken the central 
values appearing in Table IV for the A^j matrix elements for both plots. 



chosen values of the reduced matrix elements are consistent with Table II and are given by, 

A 2 J = (25, -1.3, 25, -1.0, 25, -1.8) T x 10" 4 GeV 2 (89) 

A\ = (25, -1.3, 25, -1.0, 25, -1.8) T x 10~ 4 GeV 2 (90) 

A% = (1.96, 1.96, 1.96, 1.96, 1.96, 1.96) T x 10" 3 GeV 2 (91) 

A% = (-0.035, 0.035, -0.035, 0.035) T GeV 2 . (92) 

The sum of the 27 and octet moment is given by the dashed curve in Fig. 9 which vanishes 
at the input scale of Q 2 = 5 GeV 2 . At smaller values of Q 2 , the dashed curve deviates 
from zero. We have multiplied the moment by Q 2 so that one can see more clearly that the 
spoiling of this cancellation is due to QCD evolution. The magnitude of the departure from 
this cancellation is well below the present CLAS experimental error, so an observation of 
this effect - should it be realized in nature - would again require significant improvements 
in experimental precision, at least for the leading moment. 

XII. CONCLUSIONS AND OUTLOOK 

In this work, we have provided the first complete computation of the logarithmic evolution 
of twist-four contributions to the flavor non-singlet, leading moment of the structure function 
F 2 . Our results can be employed to analyze future deep-inelastic scattering experiments 
carried out in the resonance regime, going beyond the simple ansatz given in Eq. (3) and 
utilized in Refs. [11, 18, 19]. In the present instance, we have used our calculation to carry 
out an illustrative phenomenological study the magnitude and Q 2 -dependence of the leading 
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FIG. 9: We have tuned the values of A\ and Aq to cancel against the flavor 27 moment at 
Q 2 ~ 5GeV 2 . The chosen values for the reduced matrix elements, shown in Eq.(89), are consistent 
with Table II. We have multiplied all curves by Q 2 to illustrate the breaking of this cancellation is 
due to QCD evolution. 



moment, drawing on existing lattice computations and phenomeno logical determinations of 
a subset of the operator matrix elements and plausible ansatz for the others. Despite lacking 
a complete set of non-perturbative QCD matrix element computations, we are able to draw 
a few broad conclusions from our study: 

• Theoretically one expects the overall scale of twist-4 contributions to the leading, flavor 
non-singlet moment to be small and consistent in magnitude with the results obtained 
from the CLAS analysis [11]. 

• The suppressed scale may result either from all matrix elements individually being 
small or a cancellation between various twist-4 contributions. 

• The CLAS data may provide early hints of a hierarchy among operator matrix el- 
ements, possibly indicating a larger magnitude for the flavor 8 than for the flavor 
27 

• Additional data on F 2 for the neutron as well as improvement in overall experimen- 
tal precision will be needed to further disentangle contributions from different flavor 
channels as well as to observe a cancellation scenario if it exists. 

• A full QCD prediction for the leading moment requires computation of several non- 
perturbative matrix elements. 
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To amplify on these remarks, we note that the opposite sign behavior in Fig. 7 allows for 
cancellation effects within twist-4. Such cancellations have been alluded to as the reason the 
higher twist moments are essentially independent of Q 2 in the resonance region [11], leading 
to a scaling behavior in the moment of F 2 (x,Q 2 ) known as Bloom-Gilman duality. Within 
past analyses, however, cancellations are typically found between the twist-4 and higher 
twist pieces in the twist expansion. Here we find evidence that even within twist-4, there 
is a possibility of cancellation independent of the higher twist matrix elements being small. 
However, we note that the twist-4 contribution to the moment is quite small over a wide 
range of Q 2 . It is well known that the elastic contributions to the moment begins to dominate 
in the lowest regions of Q 2 in Fig. 8, see Ref.[10]. As pointed out there, an analysis of the 
data in terms of the OPE in this region is hampered by the fact that the low n moments are 
mainly sensitive to the inclusion of elastic contributions, clouding interpretations in terms 
of higher twist effects. 

This issue can be viewed as being due to the relative sizes of the matrix elements. We 
note that the scale of both the 27 and 8 moments is set by the matrix elements of the twist- 
4 operators, and given different, possibly larger matrix elements, one may see adequate 
variation in the moment at a higher scale in Q 2 . The current experimental precision in 
extracting the leading moment for higher twist does not allow a determination of the separate 
flavor channels computed here. Thus probing cancellation effects would require higher level 
of experimental precision. From the theoretical perspective, one need not stop at a leading 
moment analysis however. 

Higher twist effects are expected to play a larger role at a given Q 2 for larger moment. 
Heuristically this is due to the fact that the resonance region weights the large- a; region more 
heavily and is thus dominated by larger-n moments. This effect is clearly illustrated in the 
CLAS data for ep scattering where the moment is decomposed in terms of the parametriza- 
tion given in Eq.(3). It is thus desirable to compute the anomalous dimension for the opera- 
tors appearing in Eq.(18) for arbitrary n since experiments are more sensitive to higher twist 
effects for larger n. The technical challenges encountered in going beyond the leading mo- 
ment are quite demanding however, both from perturbative QCD viewpoint and also from 
a lattice QCD perspective. One of the main advantages in performing a leading moment 
analysis is due to the already existing lattice estimates of the tree level 4Q operators, and 
we know of no lattice computation of these matrix elements beyond leading spin that exists 
in the literature. Due to these technical challenges, we leave a higher spin analysis of the 
type presented in here for future work. 
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Appendix A: Anomalous Dimension Matrices 



We give here the anomalous dimension matrices for operators introduced in section VI. 
When computing the one loop corrections to these operators, each operator was inserted 
into a Green's function with the same number of external legs as the operator itself, and 
dimensional regularization was employed where d = 4 — 2e. For the 2Q Feynman diagrams 
shown in Fig. 29, we have calculated each using the background field method and have chosen 
the Feynman gauge for all calculations [41]. Incorporating operator mixing, the anomalous 
dimension is, in general, a matrix in flavor space 

Following the notation in section VII the 27-plet anomalous dimension matrix is 
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The anomalous dimension matrix for the octet sector has the schematic form 
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For L/ = i, the only non-zero sector is the four-quark mixing back to four-quark piece. 
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For M/=i, we find its anomalous dimension to be 
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The anomalous dimension for 07=1,0 vanishes. For N/=i the anomalous dimension has a 
schematic form 



N 



1=1 



N 



4Q^4Q 







M 2Q^4Q T x J 2Q->2Q 
^/=1 1^/=1 / 10xl0 



(A7) 



Where the non- vanishing sectors of this matrix are, 
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The matrix L/ =0 has only one, non-vanishing piece 
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Similarly M/=o has only one sector which does not vanish, it is 
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Appendix B: Wilson Coefficients 

We list here the full expressions for the Wilson coefficients derived in [20]. Starting with 
the forward Compton amplitude, 
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In the second line above, we have parametrized the symmetric and traceless matrix elements 
in terms of the nucleon four-momenta, 



. (N|0^= 4 M JN) = (Pmx • • • P*. - traces) . 

The sum over n is even, i distinguishes the type of operator, and n denotes the spin of the 
operator, we have also used the shorthand expressions, 
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The coefficients of the tensor contribute to Fl and whereas the coefficients of d uu con 



tribute to F 2 . The full expressions for and are [22] 
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Where 



for n — 2, k — I — for example, we find e.g. 
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Using the parametrization in Eq.(53) reproduces the expression given in section VIII. 
Appendix C: Operator Reduction 

Organizing the pole structure of the one loop corrections to the 2Q operators is a non- 
trivial task due to the appearance of various non-canoical, gauge-invariant operators appear- 
ing at the one loop level. The following list of identities was used repeatedly to remove each 
non-canonical operator in favor of canonical ones. 



1. 








2. 






-i(D"j a - D a ^)il) 


3. 








4. 


D 2 






5. 






^-([D^ [Dp, D v \] + D 2 D U - D V D 2 ) 
2ig 


6. 


[A*> Dp] 






7. 









8. 






gr a ^ /7 p t>/ 



As a simple example, one such non-canonical operator encountered in the one loop analysis 
of Qf' 0) is A^A^ ^^5 v D 2 ip. Applying identity 4. and identity 1. we find this operator is 
easily expressed in terms of a linear combination of canonical operators 
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An exhaustive proof that all non-canonical, gauge invariant operators can be reduced to 
canonical form can be found in [22]. 



Appendix D: Renormalization Group Running 



The /i-dependence of the coefficient functions on Q 2 is given by the solution to the RG 
equation 
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the standard solution is given by 
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And for a one-loop analysis, the T-ordering can be dropped when evaluating the integral. 
Using the operator rescaling outlined in section V in the strong coupling, we may write 
— 9 2 (t)dji, and dropping the subscript n, r for simplicity 
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where R is a rotation matrix that diagonalizes the anomalous dimension, and we have made 
use of the following relations 



Appendix E: Matrix Elements 



In this section we list the values of all matrix elements used to predict Mj l z i \ T=4 {Q 2 ). In 
Table III and IV we include a brief description of the assumptions and denote the flavor 
structure of each operator for clarity. Following the notation at the end of section X, we list 
the reduced matrix elements (A±) the four-quark operators, 

= A J ' 2 ± 2A jA + A j ' 6 (El ) 

Ai = A jA ± 2^' 3 + A j ' 5 (E2) 

where the index j denotes the flavor representation and c denotes the presence of an SU(3) 
color generator. 

To make contact with the operators listed in Eq.(18), we have written the reduced ma- 
trix elements using the notation (<5 J ), where the index j denotes basis operator of type j 
appearing in Eq.(18). For all operators appearing in Table IV, we have set n = 2 and specify 
the k values for operators Q 7,8 . 
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TABLE III: Reduced matrix elements for the four-quark operators of each flavor representation. 
All values have been computed at the input scale 5 GeV 2 and the third column summarizes the 
methods used to compute the final value. 
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Appendix F: Conventions and Feynman Rules 

In this section, as elsewhere in the paper, we contract all free Lorentz indicies with a light 
light vector A^, where A 2 = to project out the symmetrized and traceless portion of each 
operator. Following [22] we write 

f a = ApF^ a 
f = A^ a F^ 
d = iA ■ D 

(Fl) 



where 



D, = d, + igr a Al (F2) 
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TABLE IV: All values of the reduced matrix elements used in computing M^ = 2 t =a {Q 2 ) a ^ an 
input scale of Qq = 5 GeV 2 . The methods and assumptions used to arrive at these values are 



summarized in the third column. 


(1=1,27) ibQibibQib 


Value (GeV 2 ) 


Method 


(Q 1 ) 


2.5 x 10~ 3 


Lattice [36], & Eq.(85) 


(Q 2 ) 


-1.3 x 10" 4 


Lattice [36], & Eq.(85) 


(Q 3 ) 


2.5 x 10~ 3 


Lattice [36], & Eq.(85) 




-1.0 x 10~ 3 


Lattice [36], & Eq.(85) 


(Q 5 ) 


2.5 x 10~ 3 


Lattice [36], & Eq.(85) 


(Q 6 ) 


-1.8 x 10~ 3 


Lattice [36], & Eq.(85) 


(1=1, 8,4) ipQtp i>Qi> 


Value (GeV 2 ) 


Method 


(Q 1 ) 


2.5 x 10~ 3 


27 = 8A, Assumption 


(Q 2 ) 


-1.3 x 10~ 4 


27 = 8A, Assumption 


(Q 3 ) 


2.5 x 10" 3 


27 =. 8 A, Assumption 


(Q 4 ) 


-1.0 x 10~ 3 


27 = 8 A, Assumption 


(Q 5 ) 


2.5 x 10" 3 


27 = 8A, Assumption 


(Q 6 ) 


-1.8 x 10~ 3 


27 = 8 A, Assumption 


(1=1, 8B) ^Q 2 ibij)ib 


Value (GeV 2 ) 


Method 




-0.0015 < (Q 1 ) < 0.01 


Phenomenology [37] 


(Q 2 ) 


-0.0015 < {Q 2 ) < 0.01 


Phenomenology [37] 


(Q 3 ) 


-0.0015 < (Q 3 ) < 0.01 


Phenomenology [37] 


(Q 4 ) 


-0.0015 < (Q 4 ) < 0.01 


Phenomenology [37] 


(Q 5 ) 


-0.0015 < (Q 5 ) < 0.01 


Phenomenology [37] 


(Q 6 ) 


-0.0015 < (Q 6 ) < 0.01 


Phenomenology [37] 


(i=i, 8C) ^qV 


Value (GeV 2 ) 


Method 




-0.05 < (Q 7 ' fc=0 ) < -0.02 


Phenomenology [37], & Eq.(87) 


(Q7,( fe =l) ) 


0.02 < (Q 7 ' k=0 ) < 0.05 


Phenomenology [37], & Eq.(87) 


(Q 8 '( fc =°)) 


-0.05 < (Q 7 ' k=0 ) < -0.02 


(Q s ) = {Q 7 }, Assumption 


( Q8,( fc =l) ) 


0.02 < (Q 7 ' k=0 ) < 0.05 


(Q 8 ) = (Q 7 ), Assumption 
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And the QCD vertices and propagators are 



o 
o 
o 
o 

o > 
o 
o 
o 

© 



"0"0"0"0"0"0"0"0"ff^ 



k 1 ,a,a 



gf abc \ g a M -k$- jk%) + g aP {k x - k x + -,k s x ) + g x p(k* - k? a ) 



„i 1.3 1 1.2 



i , l 



-3 ,„2 



4 



(F3) 



(F4) 



a,a p,c 



-iS. 



k 2 



(i-0 



(F5) 



+ m/) 



A; 2 — ?7ij + ie 



(F6) 



Where, in our calculations, we take rrif = 0, and choose £ = 1. The Feynman rules for the 
single gluon, 2Q operators 




q,a,fi 



' Ql =Q = 2iA a r a e a ^ 7 V q u A- Pl 
Q 7 k=1 = -2iA a r a e a ^ 7 V q u A ■ p 2 



(F7) 




q,a,fi 



Ql=o= -(A-?7^-^A m )t*A- Pi 
Qk=i= (A-q lfl -^)T a A-p 2 



(F8) 



Feynman rules for the two-gluon, 2Q operators 
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Qk=o 



Ql=i 



ig 7/375 e^ av A M {2A u T a T b ql + ir c f cah g^ v A ■ Pl } + 
i9 7/375 ^ v A„ {2A a rV g 2 , + it c f cba g au A ■ Pl } 



-i 5 7/375 e^"" A„ {2A w rV g J + ir c f cba 9uJV A ■ p 2 ) 
ig 7/375 ^ v A M {2A«T a r 6 g 2 + iT cfcab gavA . p2 } 



(F9) 




Qt=o = g {(r a r 6 A-)(A • gV - ^A a ) - (igr c f cab )(A a j u) A ■ Pl )} + 
g {{r b T a A a )(A • gV - ^ 2 A W ) - (igT c f cba )(A UJ >y a A-p 2 )} 

Ql =l = g{-( T a T b A a )(A • gV - ^ 2 A W ) + {igr c f cab )(A a ^ A ■ p 2 )} + 
5 |_( T 6 T °A a ')(A • gV - ^A a ) + (igT c f cba )(A u -f a A ■ p 2 )} 

(F10) 
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